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ABSTRACT 

Using improved doppler tomographic signal-analysis techniques we have carried out 
a deep search for starlight reflected from the giant planet orbiting the star r Bootis. 
We combined echelle spectra secured at the 4.2 ni William Herschel telescope in 1998 
and 1999 (which yielded a tentative detection of a reflected starlight component from 
the orbiting companion) with new data obtained in 2000 (which failed to confirm the 
detection). The combined dataset comprises 893 high resolution spectra with a total 
integration time of 75'i'"32""" spanning 17 nights. We establish an upper limit on the 
planet's geometric albedo p < 0.39 (at the 99.9% significance level) at the most prob- 
able orbital inclination i ~ 36°, assuming a grey albedo, a Venus- like phase function 
and a planetary radius Rp — 1.2i?j„p. We are able to rule out some combinations of 
the predicted planetary radius and atmospheric albedo models with high, reflective 
cloud decks. Although a weak candidate signal appears near to the most probable 
radial velocity amplitude, its statistical significance is insufficient for us to claim a 
detection with any confidence. 
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1 INTRODUCTION 

Two years after the discovery of a planetary companion to 
51 Pegasi by Mayor & Queloz (1995) came the identification 
of a similar object orbiting the F7V star r Bootis (Butler 
et al. 1997). Among the more surprising features of these 
discoveries were their short orbital periods, putting both 
objects very close to their parent stars. Indeed, of the 102 
extra-solar planets currently known, 19 reside within 0.1 AU 
of the parent star. In addition to posing many theoretical 
questions as to how they came to be there, the existence of 
these objects in such close orbits does open the possibility 
of their detection through the reflection of the host stellar 
spectrum. 

Shortly after the r Bootis b detection, two teams (Char- 
bonneau et al. 1999; Cameron et al. 1999) initiated a spec- 
troscopic search for the reflected light component of the or- 
biting planet. Charbonneau et al. conducted 3 nights ob- 
servation in March 1997 using the HIRES echelle spectro- 
graph mounted on the Keck I 10m telescope on Mauna 
Kea, Hawaii. Although the restricted spectral range 465.8 
- 498.7 nm provided a non detection, the signal-to-noise ra- 
tio {SNR ~ 1500) was sufficient to impose a relative re- 
flected flux limit fp/ft, < 5 X 10~^, assuming a grey albedo 
reflection of the stellar spectrum. This implies a geometric 
albedo limit p < 0.3 over the spectral range investigated, 
assuming a planetary radius of 1.2 RjupUer- Cameron et al. 



obtained 10 nights of data during 1998 and 1999, with the 
Utrecht echelle spectrograph (UES) on the 4.2 m William 
Herschel telescope. By using a least-squares deconvolution 
(LSD) technique (Donati et al. 1997) on ~ 2300 spectral 
lines in the range 385 - 611 nm, they identifled a proba- 
ble refiected-light feature with fp/ f, ^ 7.8 x 10~^ (Cameron 
et al. 1999). This detection indicated an orbital velocity am- 
plitude Kp = 74±3 km s~^ which, when combined with the 
planet's orbital velocity Vp = 152 km s^^, suggested an in- 
clination for the system i = 29° and a radius Rp — 1.8Rj, 
assuming a grey geometric albedo p = 0.55. A bootstrap 
Monte Carlo analysis gave a 5% probability that the fea- 
ture was an artefact of noise. Subsequent observations over 
7 nights in March-May 2000, however, failed to confirm the 
detection (Collier Cameron et al. 2001). 

Here we report the results of a new, deep search for 
the reflected light signal via a full re-analysis of the WHT 
data, combining all 17 nights of echelle spectra obtained 
during the 1998, 1999 and 2000 observing seasons. In Sec- 
tion 2 we use the measured system parameters to determine 
prior probabilities for the planet's orbital velocity amplitude 
and the fraction of the star's light that it intercepts. In Sec- 
tion 3 and 4 we describe the acquisition and extraction of the 
echelle spectra. Section 5 details the methods used to extract 
the planetary signal from the data, by combining the pro- 
flies of thousands of stellar absorption lines recorded in each 
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Table 1. System paraiiict(-rfs for t Bootis and its planetary eoinpanion 



Parameter Value (Uncertainty) References 



Star : 






Spectral Type 


F7V 


Fuhrmann, PfeifTer & Bernkopf (1998); Gonzalez (1998) 


Mv 


4.496 (0.008) 


Puhrmann, PfeifTer & Bernkopf (1998); Gonzalez (1998) 


Distance (pc) 


15.6 (0.17) 


Perryman, Lindegren & Kovalevsky (1997) 


TEff (K) 


6360 (80) 


Fuhrmann, Pfeiffer & Bernkopf (1998); Gonzalez (1998) 


M»(Mq) 


1.42 (0.05) 


Fuhrmann, Pfeiffer & Bernkopf (1998); Gonzalez (1998) 




1.48 (0.05) 


Puhrmann, Pfeiffer & Bernkopf (1998); Gonzalez (1998) 


[Fe/H] 


0.27 (0.08) 


Puhrmann, Pfeiffer & Bernkopf (1998); Gonzalez (1998) 


vsini (km s ^) 


14.9 (0.5) 


Henry et al. (2000) 


Age (Gyr) 


1.0 (0.6) 


Puhrmann, Pfeiffer & Bernkopf (1998) 


Planet : 






Orbital Period Porb (days) 


3.31245 (0.00003) 


Marcy, private communication 


Transit Epoch Tq (JD) 


2451653.968 (0.015) 


Marcy, private communication 


KJjns-^) 


469 (5) 


Butler et al. (1997) 


a (AU) 


0.0489 


Butler et al. (1997)(revised for this paper) 


Mp sin i(Mj„p) 


4.38 


Butler et al. (1997) (revised for this paper) 



echellogram into a time series on which we use a matched- 
filter method to measure the strength of the reflected-light 
signal. Section 7 describes how we test and calibrate the 
analysis through the addition of a simulated planetary sig- 
nal. Finally Section 8 makes theoretical assumptions about 
the size and atmospheric composition of r Bootis b in or- 
der to place upper limits on its respective geometric albedo 
and radius. We also discuss the plausibility of a candidate 
reflected-light signature that appears in the data close to 
the most probable velocity amplitude and signal strength. 



2 SYSTEM PARAMETERS 

r Bootis (HD 120136, HR 5185) is a late-F main sequence 
star with parameters as listed in Table 1. High precision ra- 
dial velocity measurements over a period of 9 years were used 
to identify a planetary companion (Butler et al. 1997) whose 
properties (as determined directly from the radial velocity 
studies or inferred using the estimated stellar parameters) 
are also summarised in Table 1. 

The equations detailed below represent a summary 
of the more comprehensive derivations described in pre- 
vious work (Cameron et al. 1999; Charbonneau et al. 
1999; Cameron et al. 2002). 

As the planet orbits its host star, some of the staxlight 
incident upon its surface is reflected towards us, producing a 
detectable signature within the observed spectra of the star. 
This signature takes the form of faint copies of the stellar 
absorption lines, Doppler shifted due to the planet's orbital 
motion and greatly reduced in intensity (~ 10~*) due to the 
small fraction of starlight the planet intercepts and reflects 
back into space. 

With our knowledge of stellar mass and the planet's 
orbital period we can estimate the orbital velocity of the 
planet Vp. The apparent radial velocity amplitude Kp of the 
reflected light is given by: 

Kp = VpSm.i = 161 sin i km s~^ (1) 
where the orbital inclination i is, according to the usual con- 



vention, the angle between the orbital angular momentum 
vector and the line of sight. 

For all but the lowest inclinations, the orbital velocity 

amplitude of the planet is substantially greater than the 
broadened widths of the photospheric absorption lines of 
the star. Hence, lines in the reflected-light spectrum of the 
planet should be Doppler shifted well clear of their stellar 
counterparts, allowing a clean spectral separation for most 
of the orbit. 

By isolating the reflected planetary signature, we are in 
effect observing the planet /star flux ratio (e) as a function 
of orbital phase {(p) and wavelength (A). 

£(a, A) s = p{\)g{a, X)^ = eo{X)g{a, A) (2) 

The phase function g{a, A) describes the variation in the 
star-planet flux ratio with illumination phase angle a. This 
is the angle subtended at the planet by the star and the 
observer, and varies according to cos a = — sini. cos<^. 

The observations measure e{a, A) over some range of 
orbital phases (f> and hence phase angles a. However, the 
signal-to-noise ratio and orbital phase coverage of the obser- 
vations is not yet adequate to define the shape of the phase 
function. Accordingly, current practice is to adopt a specific 
phase function in order to express the results in terms of 
the planet /star flux ratio that would be seen at phase angle 
zero: 

eo{X)=p{X)^ (3) 

where p(A) is the wavelength dependent geometric albedo. 

Since a is tightly constrained by Kepler's third law to 
a = 0.0489 (M,/1.42Af0) AU, the measurements of eo(A) 
measure the product Rp ■y/p(A). 

While the phase function of a Lambert sphere might be 
the simplest form to adopt for g{a,X), we prefer to adopt 
a phase function that resembles those for the cloud-covered 
surfaces of planets in our own solar system. Jupiter and 
Venus appear to have phase functions that are more strongly 
back-scattering than a Lambert sphere. Photometric studies 
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of Jupiter at large phase angles from the Pioneer and subse- 
quent missions have shown (Hovenier 1989) that the phase 
function for Jupiter is very similar to that of Venus, despite 
their very different cloud compositions. As a plausible alter- 
native to the Lambcrt-sphoro formulation, we use a poly- 
nomial approximation to the empirically determined phase 
function for Venus (Hilton 1992). The phase-dependent cor- 
rection to the planet's visual magnitude is approximated by: 

Am(a) = 0.09(q;/100°) -|- 2.39(q;/100°)^ - 0.65(a/100°)^(4) 



Prior probability density map 



so that 

g{a) = 10 ^ ' 



(5) 



In the event of any planetary detection, a careful anal- 
ysis of the data should allow us to determine the following 
information: 

(i) Kp, the planet's projected orbital velocity, from 
which we obtain the orbital inclination of the system and 
hence the mass of the planet, since Mp sini is known from 
the star's Doppler wobble. 

(ii) eo, the maximum flux ratio observed, with which 
we can constrain the planet's radius since eo = p(Rp/af', 
where p is the geometric albedo of the planetary atmosphere. 
Alternatively we can adopt a theoretical radius to constrain 
the albedo of a given atmospheric model. 



2.1 Rotational broadening 

The rotational broadening of the direct starlight and chro- 
mospheric Ca II H & K emission flux suggest that the star's 
rotation is synchronised with the orbit of the planet (Bal- 
iunas et al. 1997; Henry et al. 2000). In a tidally locked 
system there is no relative motion between the surface of 
the planet and the surface of the star, so the planet will 
reflect a non-rotationally broadened stellar spectrum, with 
typical line widths dominated by turbulent velocity fields in 
the stellax photosphere. These motions were estimated by 
Baliunas et al. (1997) to be of the order ~ 4 km s~^. Any 
absorption lines attributed to the planet's atmosphere are 
thus likely to be much narrower than the stellar lines. 



2.2 Orbital Inclination 

In the first instance we can rule out inclinations i > 80° 

due to the absence of transits in high-precision photometry 
(Henry et al. 2000). Furthermore, if we assume the star's 
rotation to be tidally locked to the planet's orbit, we can 
use the projected equatorial rotation speed of the host star 
vsini = 14.9 ± 0.5 km s"^ (Henry et al. 2000) to loosely 
constrain the orbital inclination to i ~ 40°. We would thus 
expect a projected orbital velocity amplitude close to Kp ~ 
100 km s-^ 



2.3 Planet Radius 

The HD 209458b transit detection of Charbonncau et al. 
(2000) provided the first confirmation of the gas giant na- 
ture of close-in extra solar planets, and yielded a radius in 
good agreement with the predictions of past and current 
interior structure models (Guillot et al. 1996; Seager & Sas- 
selov 1998; Marley, Gelino & Stephens 1999; Burrows et al. 
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Figure 1. The greyscale shows the prior joint probability density 
function (PDF) for projected orbital velocity Kp and the squared 
ratio (Rp/a)^ of the planet radius to the orbit radius, based on the 

measured system parameters. Darker shades in the greyscale de- 
note greater probabilities of the planet having the corresponding 
combination of {Rp/a)'^ and Kp. The PDF shows that the planet 
is most likely to have Kp ~ 94 km s~^ and (Rp/a)^ ~ 1.44x 10~*. 



2000; Seager, Whitney & Sasselov 2000). In short, the planet 
radius evolves with time and depends on the planet mass. 
For our purposes, this defines a range of theoretically plau- 
sible radii at each possible value of the planet mass. The 
range of possible planet radii was computed specifically for 
r Bootis b by Burrows et al. (2000), allowing for uncertain- 
ties in the orbital inclination and hence the planet's mass. 
Their radiative-eonvective gas giant models predict upper 
limits on the planet's radius of 1.58 Rjup for Mp = 7 Mjup, 
and 1.48 Rjup for Mp = 10 Mjup. 

2.4 Prior Estimates of System Parameters 

In searching for a faint reflected-light signature from a planet 
with an unknown orbital inclination, it is useful to know 
in advance how the planet's observable properties ought to 
depend on the orbital inclination. We do this by construct- 
ing the a priori probability density functions for the vari- 
ous observable properties of the planet described in Table 1. 
This helps us to determine whether any faint candidate re- 
flection signature is physically plausible, given our existing 
knowledge of the system's parameters. We do not want to be 
guided too closely by theory, but values outside the plausible 
ranges would pose difficulties for current thinking. 

Fig. 1 shows the probability distributions for the 
planet's radial velocity amplitude Kp and the quantity 
{Rp/af , based on a Monte Carlo simulation using the ex- 
pressions given in Section 2. We assume Gaussian distri- 



© 0000 RAS, MNRAS 000, 000-000 



4 C.Leigh, A.Collier Cameron, K.Horne, A.Penny and D.James 



Table 2. Journal of observations. The UTC mid-times and orbital phases are shown for the first and last spectral exposures secured on 
each night of observation. The number of exposures is given in the final column. 



UTC start Phase UTC End Phase Exposures 



1 onQ 


Apr 


no 

uy 




U 


yl or; 
4zO 


1 noQ A 1 n 
iyycs Apr iu 


Uo:o ( :Uo 


U 


Ki n 
Oiy 


1 n7 


1 OQQ 


Apr 


1 n 
iU 


oo.n/i ./in 


U 


1 ZD 


1 QOQ A T-vv 1 1 

lyyo Apr ii 


UD.zU.04 


U 


ooU 


1 1 Q 


1 QQQ 

iyy© 


Apr 


1 1 


oo.no.oQ 


U 


noQ 
Uzy 


1 QQQ A r\T 1 O 

lyycs Apr iz 


U0.0c5.44 


U 


1 07 
iz I 


O-L 


1 OQQ 

lyyo 


Apr 


1 Q 


OQ-no- C^yl 


U 




1 QOQ A T-iT- 1 A 

lyyo Apr 14 


UO.o ( .UO 


n 
U 


( oo 


OO 


1 ooo 
lyyy 


Apr 


no 


oo-nfi-/i 

ZZ.UD.40 


n 
U 


/I OQ 

4yo 


lyyy Apr uo 


UD.Uo.4 i 


U 


oy4 


ZO 


1999 


Apr 


25 


21:43:17 


Q 


441 


1 QQQ Ar»r 9fi 




Q 


539 


40 


1999 


May 


05 


21:56:59 





463 


1999 May 06 


04:51:57 





550 


60 


1999 


May 


25 


20:59:45 





488 


1999 May 26 


03:56:18 





576 


51 


1999 


May 


28 


20:52:23 





390 


1999 May 29 


03:04:35 





467 


47 


1999 


Jun 


04 


20:22:34 





497 


1999 Jun 05 


00:07:23 





548 


23 


2000 


Mar 


14 


23:14:49 





271 


2000 Mar 15 


06:54:52 





366 


48 


2000 


Mar 


15 


22:51:55 





567 


2000 Mar 16 


06:49:04 





669 


45 


2000 


Mar 


24 


22:21:55 





278 


2000 Mar 25 


06:47:50 





385 


34 


2000 


Apr 


23 


20:38:59 





311 


2000 Apr 24 


05:07:31 





420 


45 


2000 


Apr 


24 


20:58:32 





619 


2000 Apr 25 


05:14:11 





724 


57 


2000 


May 


13 


20:47:01 





349 


2000 May 14 


03:12:08 





434 


41 


2000 


May 


17 


20:24:14 





557 


2000 May 18 


03:52:37 





649 


43 



butions for the measured stellar mass (1.42 ± 0.05Mq), rar 
dius (1.48 ± 0.05/?©), stellar reflex velocity (469 ± 5 m s"^) 
and the radius of the planet. The planet radius and un- 
certainty range were obtained from the theoretical mass- 
radius relations (Guillot et al. 1996; Burrows et al. 2000) 
described above. Both works find the most probable radius 
to be 1.2 ± Q.lRjup, assuming r Bootis b to have an age of 
1 Gyr. 

Furthermore, we Eissume the star's rotation is tidally 
locked to the planet's orbit. We thus generate a distribution 
of sini values based on Gaussian distributions for the pro- 
jected stellar rotation speed {v sini — 14.9±0.5 km s^^) and 
a stellar rotation period (3.3 ±0.1 days) closely bound to 
the orbital period of the planet. A further restriction applies 
where the tidal synchronisation timescale for the primary's 
rotation, 

r.,„e^l.2(^) {-^y years (6) 

is longer than the main-sequence liftime of the host star 
Tms — W^'^{Mt/Mr.))^'^ years, in which case we reject that 
model from the Monte Carlo analysis. 

The resulting probability map (Fig. 1) shows a Kp dis- 
tribution centred on ~ 94 km s~^ , with the most likely value 
of {Rp/af ~ 1.44 X 10"". The projection of this PDF on 
to the orbital velocity axis defines the region of parameter 
space in which we can be confident that a detection would 
occur if it were present in the data, given our prior knowl- 
edge of the system parameters. We use the Kp projection 
of the PDF in the subsequent analysis to test the plausi- 
bility of any candidate features which appear in the data, 
by modifying the posterior probability distribution to assess 
the false alarm probability (see Section 8.2.1). Secondly, our 
upper limits on co arc sensitive to the orbital inclination, 
so we adopt the most probable Kp in order to determine 
the most plausible upper limits on the planet's radius and 
albedo. 



For any given albedo model we can also use the data 
to determine upper limits on {Rp/af instead of the opposi- 
tion flux ratio eo. The projection of the PDF onto (Rp/af 
thus allows us to compare the effective reflection area of the 
planet directly with model predictions. Unlike the projec- 
tion on to Kp, however, the prior probability distribution 
for {Rp/af plays no role in assessing the plausibility or oth- 
erwise of a candidate detection. 



3 OBSERVATIONS 

We observed r Bootis during 1998, 1999 and 2000 using the 
Utrecht Echelle Spectrograph on the 4.2 m William Herschel 
Telescope at the Roque de los Muchachos Observatory on La 
Palma. The detector was a single SITe 1 CCD array contain- 
ing some 2048 x 2048 13.5-/im pixels. The CCD was centred 
at 459.6 nm in order 124 of the 31 g mm~^ echelle grating, 
giving complete wavelength coverage from 407.4 nm to 649.0 
nm with minimal vignetting. The average pixel spacing was 
close to 3.0 km s~^, and the full width at half maximum 
intensity of the thorium-argon arc calibration spectra was 
3.5 pixels, giving an effective resolving power R = 53000. 

Table 2 lists the journal of observations for the 17 nights 
of data which contribute to the analysis presented in this pa- 
per. In the first year (1998) the stellar spectra were exposed 
between 100 and 200 seconds. For 1999 and 2000, the stellar 
spectra were exposed for between 300 and 500 seconds, de- 
pending upon seeing, in order to expose the CCD to a peak 
count of 40000 ADU per pixel in the brightest parts of the 
image. A 450-s exposure yielded about 1.2 x 10® electrons 
per pixel step in wavelength in the brightest orders in typical 
(1 arcsec) seeing after extraction. We achieved this with the 
help of an autoguidcr procedure, which improves efficiency 
in good seeing by trailing the stellar image up and down 
the slit by ±2 arcsec during the exposure to accumulate the 
maximum S:N per frame attainable without risk of saturar 
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tion. Note that the 450-s exposure time compares favourably 
with the 53-s readout time for the SITe 1 CCD in terms of 
observing efficiency - the fraction of the time spent collect- 
ing photons is above 90%. Following extraction, the S:N in 
the continuum of the brightest orders is typically 1000 per 
pixel. 



4 SPECTRUM EXTRACTION 

One-dimensional spectra were extracted from the CCD 
frames using an automated pipeline reduction system built 
around the Starlink ECHOMOP and FIGARO packages. 
Nightly flat-field frames were summed from 50 to 100 frames 
taken at the start and end of each night, using an algo- 
rithm that identified and rejected cosmic rays and other noii- 
repeatable defects by comparing successive frames. Due to 
physical movement of the chip mounting between and dur- 
ing observation runs, it was found that the level of noise was 
reduced by the use of nightly flat flelds rather than master 
flat fields for the entire year's observations. 

The initial tracing of the echelle orders on the CCD 
frames was performed manually on the spectrum of r Boo- 
tis itself, using exposures taken for this purpose without 
dithering the star up and down the slit. The automated ex- 
traction procedure then subtracted the bias from each frame, 
cropped the frame, determined the form and location of the 
stellar profile on each image relative to the trace, subtracted 
a linear fit to the scattered-light background across the spa- 
tial profile, and performed an optimal (profile and inverse 
variance-weighted) extraction of the orders across the full 
spatial extent of the object-plus-sky region. Nightly flat- 
fleld balance factors were applied in the process using the 
50 to 100 frames obtained at the start and end of each night 
of observations. In all, 55 orders ( orders 88 to 142 ) were 
extracted from each exposure, giving full spectral coverage 
from 407.4 to 649.1 nm with good overlap. 



5 EXTRACTING THE PLANET SIGNAL 

For a bright, cloudy model planet with p = 0.4 and Rp = 
1.2Rjupi we expect the flux of starlight scattered from the 
planet to be no more than one part in 18000 of the flux 
received directly from r Bootis itself, even at opposition 
(a = 0). In order to detect the planet signal, we flrst subtract 
the direct stellar component from the observed spectrum, 
leaving the planet signal embedded in the residual noise 
pattern. A detailed description of this procedure is given in 
Cameron ct al. (2002) Appendix A. The planet signal con- 
sists of faint Doppler-shifted copies of each of the stellar ab- 
sorption lines. After cleaning up any correlated fixed-pattern 
noise remaining in the difference spectra (sec Cameron et al. 
2002 Appendix B), we then create a composite residual line 
profile, by fitting to the thousands of lines recorded in each 
echellogram (Cameron et al. 2002 Appendix C). Finally we 
use a matched-filter analysis (Cameron et al. 2002 Appendix 
D) to search for features in the time-series of composite 
residual profiles whose temporal variations in brightness and 
radial velocity resemble those of the expected rcflccted-light 
signature. For an assumed albedo spectrum p(A) and orbital 



velocity amplitude Kp, the fit of the matched filter to the 
data measures {Rp/af . 



6 ANALYSIS CHANGES 

In this new analysis of the r Bootis data, we have made the 

following significant changes to the processing undertaken 
for the original Cameron ct al. (1999) paper - 

(i) The inclusion of the year 2000 data, which adds seven 
nights' data taken at optimally-illuminated orbital phases to 
the analysis. 

(ii) Full re-extraction of all three years' data, again us- 
ing optimal methods, and providing an increase in the spec- 
tral range by two echelle orders or ~ 15 nm. 

(iii) The use of nightly flat-field frames in the extraction 
routine, rather than the previous whole year flat-fields. Post 
extraction analysis showed a ^ 4% reduction in noise. 

(iv) Increases in computational processing power over 
the intervening two years has allowed the analysis to be con- 
ducted on individual echelle frames, rather than having to 
co-add the spectra into groups of four prior to the deconvo- 
lution and matched-filter analysis. 

(v) The inclusion of a Principal Component Analysis 
routine (PCA), as detailed in Cameron et al. (2002) Ap- 
pendix B, to remove correlated fixed-pattern noise that was 
appearing in the difference spectra (i.e. raw spectra - stellar 
template frames). 

(vi) The use of a more stringent calibration technique, 
described at Appendix A, to quantify and correct for the 
fraction of the planetary signal lost during the stellar sub- 
traction, deconvolution and PCA routines. With this we pro- 
duce shallower but more realistic upper limits than were 
stated by Collier Cameron et al. (2001), who assumed no 
loss of signal. 



7 SIMULATED PLANET SIGNATURES 

We verified that a faint planetary signal is preserved through 
the above sequence of operations in the presence of realistic 
noise levels, by adding a simulated planetary signal to the 
observed spectra. We also use the simulated signal to cali- 
brate the strength of any detected signal (Appendix A) . The 
simulations were based on the assumption that the planet's 
rotation is close to being tidally locked, always keeping the 
same face towards the star. The resulting broadening of the 
spectral lines is therefore dominated by convective motions 
on the star's surface, estimated at ~ 4 km s~^ (Baliunas 
et al. 1997). For our sirrmlations we chose to use the slowly 
rotating giant star HR 5694, observed on several nights in 
1999. HR 5694 is a F7III spectral type of similar temper- 
ature and elemental abundance to r Bootis, but with an 
estimated vsvai ~ 6.4 ± 1 km s"^, making it well suited to 
represent the refiected starlight (Baliunas et al. 1997). 

For any assumed axial inclination, the phase angle and 
line-of-sight velocity are known at all times. The simulation 
procedure simply consists of shifting and scaling the spec- 
trum of HR 5694 according to the orbit and phase function, 
co-multiplying it by an appropriate geometric albedo spec- 
trum, and adding it to the observed data. To ensure a strong 
signal we used a simulated planet of radius 1.4 iijup and 
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Phased dynamic spectrum Phased dynamic spectrum 




Velocity km/s 

Figure 2. Time series of deconvolved profiles derived from the 
original WHT spectra, and secured over 17 nights observations in 
1998, 1999 and 2000, but with the addition of a simulated planet 
signal at an inclination of 60°. The injected signal is that of a 
planet with geometric albedo p = 0.5 and radius lARj^p. The 
planetary signature appears as a dark sinusoidal feature crossing 
from right to left as phase increases and centred on the superior 
conjunction at phase 0.5. 




-200 -100 100 200 

Velocity km/s 

Figure 3. Time series of deconvolved profiles derived from the 
original WHT spectra, and secured over 17 nights observations 
in 1998, 1999 and 2000. The observations assume a grey albedo 
spectrum without the addition of a synthetic planet signal. The 
greyscale runs from black at — 10~* times the mean stellar contin- 
uum level, to white at +10~*. The velocity scale is in the reference 
frame of the star. 



wavelength-independent geometric albedo p = 0.5, which 
when viewed at zero phase angle should give a planet-to- 
star flux ratio eo — fp/f* = 0.98 x lO""*. We have chosen a 
planetary radius greater than that expected by theory so as 
to provide a simulated input signal strong enough to return 
an unambiguous detection. 

The resulting time series of deconvolved line profiles, 
shown in Fig. 2, demonstrates how the simulated planet sig- 
nal is recovered after the extraction process, with the plane- 
tary signal clearly visible as a dark sinusoidal feature cross- 
ing from right to left between phases 0.25 and 0.75. The 
weakening of the simulated planetary signature near quadra- 
ture is caused mainly by the phase function. The signal is 
further attenuated near quadrature by the way in which the 
templates are computed: since the planet signature is nearly 
stationary in this part of the orbit, some of the signal will be 



removed along with the stellar profile if many observations 
are made in this part of the orbit. 

Figs. 2 and 3 both show a "barber's-pole" pattern of dis- 
tortions in the residual stellar profiles at low velocities. The 
phase variation in these undulations appears consistent with 
sub-pixel shifts in the position of the spectra with respect to 
the detector over the course of the night. Fortunately they 
only affect a range of velocities at which the planet signature 
would in any case be indistinguishable from that of the stax. 

The relative probabilities of the fits to the data for dif- 
ferent values of the free parameters Rp/a and Kp are given 
by 

P{Kp,Rp/a)cxeKp{-xV2), (7) 
where 
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Kp l<m/sec 

Figure 4. Relative probability map of model parameters Kp and 
log(eo) = logp(iip/a)^ for a simulated planet signature with grey 
albedo p = 0.5, Rp = lARj^p and orbital inclination of 60 °. 

The contours show the confidence levels at which candidate detec- 
tions can be ruled out as being caused by spurious alignments of 
non-Gaussian noise features. From top to bottom, they sfiow tfic 
99.9%, 99.0%, 95.4% and 68.4% confidence limits. The synthetic 
planet signature is detected well above the 99.9% confidence limit. 



Kp l<m/sec 

Figure 5. Relative probability map of model parameters Kp and 
log(eo) = logp(Rp/a)2, derived from the WHT/UES observa- 
tions of r Bootis, assuming a grey albedo spectrum. The greyscale 

denotes the probability relative to the best-fit model, increasing 
from for white to 1 for black. A broad candidate feature appears 
close to the 99.0% confidence contour, with projected orbital ve- 
locity amplitude Kp = 97 km s~^. 



This is conveniently displayed in greyscale form as a function 
of Kp and log(eo) = logp(i?p/a)^. In Fig. 4 we show the 
map for the simulated observations, with the probabiUties 
normalised to the most probable value in the map. 

The signal of the synthetic planet appears as a compact, 
dark feature at Kp = 139 km and log(eo) = —4.01, 
i.e. eo ~ 0.98 x 10~*. This most probable combination of 
orbital velocity and planet radius yields an improvement 
Ax^ = 84.6 with respect to the value obtained assuming no 
planet is present (Fig. 6). 

To set an upper limit on the strength of the planet sig- 
nal, or to assess the likelihood that a candidate detection is 
spurious, we need to compute the probability of obtaining 
such an improvement in by chance alone. In principle 
this could be done using the distribution for 2 degrees 
of freedom. In practice, however, the distribution of pixel 
values in the deconvolved difference profiles has extended 
non-Gaussian tails that demand a more cautious approach. 

Rather than relying solely on formal variances derived 
from photon statistics, we use a "bootstrap" procedure to 
construct empirical distributions for confidence testing, us- 
ing the data themselves. In each of 3000 trials, we random- 
ize the order in which the 17 nights of observations were 
secured, then we randomise the order in which the obser- 



vations were secured within each night. The re-ordered ob- 
servations are then associated with the original sequence of 
dates and times. This ensures that any contiguous blocks 
of spectra containing similar systematic errors remain to- 
gether, but appear at a new phase. Any genuine planet sig- 
nal present in the data is, however, completely scrambled 
in phase. The re-ordered data are therefore as capable as 
the original data of producing spurious detections through 
chance alignments of blocks of systematic errors along a sin- 
gle sinusoidal path through the data. We record the least- 
squares estimates of log(eo) ~ logp(i?p/a)^ and the associ- 
ated values of as functions of Kp in each trial. 

The percentage points of the resulting bootstrap distri- 
bution are shown as contours in Figs. 4 and 5. Prom bot- 
tom to top, these contours give the 68.4% 95.4% 99.0% and 
99.9% bootstrap upper limits on the strength of the planet 
signal. The 99.9% contour, for example, represents the value 
of log(eo) that was only exceeded in 3 of the 3000 trials at 
each Kp. 



8 RESULTS AND DISCUSSION 

The results of this analysis appear on the relative probabil- 
ity map of model parameters Kp and log(eo) ~ \ogp{Rp/a)^ , 
shown at Fig. 5. The calibrated confidence levels allow us 
to achieve our primary aim of constraining the radius and 
albedo of the planet. However, there exists a significant can- 
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Figure 6. The upper panel shows the optimal scaling factor 

(Rp/a)-^ plotted against orbital velocity amplitude Kp, assuming 
a grey albedo spectrum for the simulated planet data. The lower 
panel shows the associated reduction Ax^ = 84.6, measured rel- 
ative to the fit obtained in the absence of any planet signal i.e. 
for (Rp/a)^ = 0. Note that only positive values of {Rp/a)^ are 
physically plausible. 

didate feature close to the 99% level that requires further in- 
vestigation. In the subsequent discussion, we therefore also 
explore the possibility that this feature could represent a 
genuine planetary detection. 

If the feature were genuine, the projected orbital ve- 
locity amplitude Kp ~ 97 (±10) km yields an orbital 
inclination of 37 (±5)°. This would be consistent with the 
star's rotation being tidally locked to the planet's orbit and 
implies a mass for r Bootis b of Mp = 7.28(±0.83)Mj„p. 

We emphasise that, although the feature appears very 
close to the peaJj of the prior probability distribution pro- 
jected onto Kp, shown in Fig. 1, there remains a distinct 
possibility that the candidate detection is a consequence of 
spurious noise and as such we should proceed with caution. 



8.1 Upper Limits on Grey Albedo 

The grey albedo model assumes that at all times the planet- 
star flux ratio is independent of wavelength. For an assumed 
planetary radius we can thus use Equation 3 to constrain 
the geometric albedo. Table 3 lists the upper limits on the 
albedo at various levels of significance, for the planetary ra^- 
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Figure 7. As for Fig. 6 but assuming a grey albedo spectrum 
for the original data without a simulated planet signal. The lower 
panel shows the associated reduction in Ax^ = 8.324, mcasiu'cd 
relative to the fit obtained in the absence of any planet signal, l^'he 
weak improvement in Ax^ at Kp = 60 km s~^ corresponds to a 
negative value of {Rp/a)^ and is therefore physically implausible. 



Table 3. Upper limits on the grey albedo for the atmosphere of 
T Bootis b assuming a radius of 1.2 Rjup- The upper limits are 
quoted for an assumed Kp ~ 94 km s~^, at the peaJf of the prior 
probability distribution for Kp. 



False Alarm p{Rp/a)^ Upper Albedo Limit 



Probability p 

0.1 % 0.561 E-04 0.39 

1.0 % 0.403 E-04 0.28 

4.6 % 0.305 E-04 0.21 



dius Rp = 1.2Rjup predicted by current theoretical models 
(Guillot et al. 1996; Burrows et al. 2000). 

The contours in Fig. 5 produced by the bootstrap sim- 
ulation constrain the maximum reflected flux ratio at oppo- 
sition to be £() < 0.561 X lO^'* at the 99.9 % confidence level, 
assuming a projected orbital velocity Kp at the peak of the 
prior probability distribution, i.e. Kp ~ 94 km s~^. This 
would limit the geometric albedo of the planet to p < 0.39. 
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Figure 8. Geometric albedo spectra for the Class V (solid line), 
isolated (dashed) and irradiated (dot-dashed) Class IV models of 
Sudarsky, Burrows & Pinto (2000), plotted over the wavelength 
range we observed. Together with a grey model of albedo p = 0.3, 
the Class V and isolated Class IV models were used to probe the 
wavelength dependence of candidate reflected-light signals. The 
geometric albedo spectrum of Jupiter (short dashes) is shown for 
comparison. 



Table 4. Upper limits on planet radius for various albedo models. 
The limits are quoted for an assumed Kp ~ 94 km s~^, at the 
pealj of the prior probability distribution for Kp . Note that the 
results for the grey albedo model are given for a geometric albedo 
of p = 0.3. Radii for other grey model albedos can be obtained 
by dividing Column 4 by -^p/O.S . 



Albedo Model 


False Alarm 


(Rp/ar 


Rp/ Rjup 




Probability 




Upper Limit 


Grey 


0.1% 


1.87 E-04 


1.37 


(P = 0.3) 


1.0% 


1.34 E-04 


1.16 




4.6% 


1.02 E-04 


1.01 


Class V 


0.1% 


1.16 E-04 


1.08 




1.0% 


0.91 E-04 


0.95 




4.6% 


0.63 E-04 


0.79 


Class IV 


0.1% 


1.50 E-04 


1.22 


(Isolated) 


1.0% 


1.13 E-04 


1.06 




4.6% 


0.87 E-04 


0.93 



We note that this is a similar result to that obtained by 
Charbonneau et al. (1999) at the same inclination. Both 

studios assume a grey albedo, Rp = 1.2Rjup, synchronous 
rotation of the star and hence a reflected version of the stel- 
lar spectrum with no rotational broadening. The candidate 
feature that appears in Fig. 5 would, if genuine, yield a grey 
geometric albedo of p = 0.32 (±0.13) for a planet of this 
radius. 

8.2 Upper Limits on Radius 

Here we investigate how atmospheric albedo models can be 
incorporated into the signal analysis to place upper limits 

on the size of the planet. A non-grey albedo model is built 
into the formation of the least-squares deconvolved profile, 
by scaling the strengths of the lines in the deconvolution 
mask by a factor equal to the geometric albedo at each line's 
wavelength. The scale factor produced by the matched-filter 
analysis is then directly proportional to (Rp/a)^ (sec Equa- 
tion 3) and is calibrated by injecting the signature of a planet 
of known radius with the specified albedo spectrum into the 
data. The method is described in detail by Cameron et al. 
(2002). 

The theoretical models we consider here arc those pro- 
posed by Sudarsky, Burrows & Pinto (2000) for a range of 
extrasolar giant planets, as shown at Fig. 8. These mod- 
els axe grouped primarily by their mass and orbital distance 
from the host star, factors which in turn influence their effec- 
tive surface temperature, surface gravity and hence radius. 
We recognise, however, that recent observations of the at- 
mosphere of HD 209458b (Charbonneau et al. 2002) suggest 
there may be less sodium absorption than predicted in these 
and other models (Brown et al. 2001; Hubbard et al. 2001). 

8.2.1 Grey Albedo Model 

At the most probable values in the prior distribution {Kp ~ 
94 km s"^), assuming a grey albedo model of p = 0.3, the 



0.1%, 1.0% and 4.6% upper limits on the planet/star flux 
ratio eo correspond to upper limits on the planet radius, 
as detailed in Table 4. With a higher assumed geometric 
albedo, the planet's radius is more strongly constrained. 

Our potential planet signal yields an improvement 
Ax^ = 8.324 over the model flt obtained assuming no planet 
signal is present (Fig. 7). We used the bootstrap simulations 
to determine the probability that a spurious feature with 
Ax^ > 8.324 could be produced by a chance alignment of 
noise features in the absence of a genuine planet signal. It 
is important to note that the bootstrap contours only give 
the false-alarm probability if the value of Kp is known in ad- 
vance, which is not the case here. The true false-alarm prob- 
ability is greater, being the fraction of bootstrap trials where 
spurious peaks at any plausible value of Kp can exceed the 
Ax^ of the candidate. If we assume that all values of Kp are 
equally likely in the range 50 km s"^ <Kp < 162 km s"\ 
the false-alarm probability is found to be 14.7% via the 
method described more comprehensively in Cameron et al. 
(2002) Appendix E. 

In practice, however, we are more likely to believe that 
a feature detected near the peak of the prior probability dis- 
tribution for Kp is genuine, than if the feature appeared at 
a velocity that was physically implausible given our existing 
knowledge of the system parameters. We can therefore use 
our prior estimation of Kp to weight the false-alarm prob- 
abilty, in the manner discussed in Cameron et al. (2002) 
Appendix E. We find from Table 5 that the false-alarm 
probability drops to 3.6% when prior knowledge of Kp is 
accounted for. For comparison, we find that a matched fil- 
ter analysis of the simulated planet data (Fig. 6) sees an 
improvement of Ax^ ~ 84.6 above the value obtained as- 
suming no planet signal is present. This is far greater than 
the Ax^ = 42.3 produced at any Kp in the bootstrap tri- 
als. The false-alarm probability is therefore substantially less 
than one part in 3000, and as such the "detection" of the 
simulated signal is secure. 
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Table 5. Projected orbital velocity peak, planet radius, Ax'^ and false-alarm probabilities (FAP) for the candidate feature, on the basis 
that it represents a genuine detection.The second FAP weights Kp in proportion to the prior probability density distribution (Fig. 1). 



Albedo Model 








FAP 


FAP 




{km s~^) 






(Uniform Weight) 


(Kp Prior) 


Grey (p = 0.3) 


97 (±8) 


1.24 ± 0.25 


8.324 


0.147 


0.036 


Class V 


95 (±8) 


1.08 ±0.19 


12.06 


0.032 


0.003 


Class IV 


90 (±8) 


1.18 ±0.20 


9.227 


0.092 


0.032 



Upper limits 

T I I I I I I I I I I I I r 
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K km/sec 

Figure 9. Relative probability map of model parameters Kp and 
log(i?p/a)^, derived from the WHT/UES observations of r Boo- 
tis, assuming the albedo spectrum to be that of a Class V roaster. 
The greyscale and contours are defined as in Fig. 4. 

8.2.2 Class V Model 

The "Class V roaster" is the most highly reflective of the 
models published by Sudarsky, Burrows & Pinto (2000). It 
is characteristic of planets with T^ff > 1500 K and /or sur- 
face gravities lower than ~ 10 m s~^, and as such is as- 
sociated with lower mass planets, such as v And b. The 
model predicts a silicate cloud deck located high enough in 
the atmosphere that the overlying column density of gaseous 
alkali metals is low, allowing a substantial fraction of incom- 
ing photons at most optical wavelengths to be scattered back 
into space. There remains, however, a substantial absorption 
feature around the Na I D lines, as shown in Fig. 8. 

We carried out the deconvolution using the same line 
list as for the grey model, but with the line strengths at- 
tenuated using the Class V albedo spectrum (see Cameron 
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Figure 10. As for Fig. 6 but assuming a Class V spectrum for the 
original data without a simulated planet signal.The lower panel 

shows the associated reduction in Ax^ of 12.06, measured relative 
to the fit obtained in the absence of any planet signal. 

et al. (2002) Appendix C). We calibrated the signal strength 

as described in Cameron ct al. (2002) Appendix D, by inject- 
ing an artificial planet signature consisting of the spectrum 
of HR 5694, attenuated by the Class V albedo spectrum 
and scaled to the signal strength expected for a planet with 
Rp = l.^Rjup. 

The form of the Class V probability map, as shown in 
Fig. 9 is similar to that encountered for the grey albedo 
spectrum. The resulting upper limits on the planet radius 
are detailed at Table 4, with the corresponding false-alarm 
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Figure 11. Relative probability map of model parameters Kp 
and log(ilp/o)^, derived from the WHT/UES observations of r 
Bootis, assuming the albedo spectrum to be that of an "isolated" 

Class IV gas giant. The candidate feature, if genuine, corresponds 
to the detection of a 1.18(±0.20) Rj planet. The greyscale and 
contours are defined as in Fig. 4. 

probabilities listed in Table 5. We find the feature produces 
a local probability maximum near Kp = 95 km s~^, with an 
improvement in Ax^ over the grey albedo model of 12.06, 
as plotted in Fig. 10. This improvement translates to a 
reduced FAP (unweighted) of 3.2%, however, the position 
of the best-fitting Kp matches the prior probability maxi- 
mum of Kp = 94 km s~^ so closely that the overall FAP 
is 0.3%, substantially lower than that obtained for the grey 
albedo case. This suggests strongly that the features in the 
data that give rise to this signal originate predominantly 
at blue wavelengths. If the candidate feature we observe 
were genuine, it would indicate a Class V planet of radius 
1.08(±0.19) Rj, which is in line with with current theory 
(Guillot et al. 1996; Burrows et al. 2000). 

8.2.3 Isolated Class IV model 

The "Class IV" models of Sudarsky, Burrows & Pinto (2000) 
have a more deeply-buried cloud deck than the Class V mod- 
els and are probably more closely applicable to r Boo b given 
its relatively high surface gravity. The resonance linos of Na 
I and K I are strongly saturated, with broad damping wings 
due to collisions with H2 extending over much of the optical 
spectrum (Fig. 8). 

We used the procedures described above to deconvolve 
and back-project the data assuming an "isolated" Class IV 
spectrum. Although this model does not take full account 
of the effects of irradiation of the atmospheric temperature- 
pressure structure, it is a useful compromise between the 
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measured relative to the fit obtained in the absence of any planet 
signal. 



Class V models and the very low albedos found with irra- 
diated Class IV models. The resulting time-series of decon- 
volved spectra is noisier than the Class V and grey-albedo 
versions, because lines redward of 500 nm contribute little 
to the deconvolution. 

The probability map (Fig.ll) derived from the boot- 
strap matched-filter analysis again shows a marginally sig- 
nificant candidate reflected-light feature, but for this albedo 
model the peak of the distribution is shifted to Kp = 
90 km s~^. The fit to the data is slightly better than in 
the grey albedo case, giving an improvement of Ax^ ~ 
9.227 over the no-planet hypothesis. The bootstrap analy- 
sis returns an unweighted false-alarm probability of 9.2%, 
but the displacement of the most probable value of Kp 
(90 km s"^) away from the prior probability maximum at 
Kp = 94 km s~^ means the overall Bayesian FAP is 3.2%, 
slightly lower than in the grey albedo case. 

The increase in noise associated with the extraction of 
the Class IV simulation produces more loosely constrained 
upper limits on the planetary radius, as set out in Table 4. 
We note that if the candidate feature were genuine, it would 
indicate a Class IV planet of radius 1.18(±0.20) Rj, still in 
line with current expectations (Guillot et al. 1996; Burrows 
et al. 2000). 
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9 CONCLUSION 

We have re-analysed the WHT echelle spectra obtained for 
the F7V star r Bootis during 1998, 1999 and 2000. By as- 
suming that 

(a) the rotation of t Bootis is tidally locked to the orbit 
of the planetary companion, suggesting a orbital inclination 
of i ~40°, and 

(b) the planet radius Rp ~ 1.2Mjup, in accordance with 
the general theoretical predictions of Guillot et al. (1996); 
Burrows et al. (2000) and with observations of the planetary 
transits across HD 209458 (Charbonneau et al. 2000), 

we are able to rule out a reflective planet with a grey albedo 
greater than p = 0.39 to the 99.9% confidence level. The 
alternative approach of adopting the specific grey (p = 0.3), 
Class V and Class IV albedo models of Sudarsky, Burrows 
& Pinto (2000) places model-dependent upper limits on the 
planetary radius. The results indicate with 99.9% confidence 
that the upper limits are 1.37 Rjup, 1-08 Rjup and 1.22 Rjup 
respectively for the three albedo models considered. 

Our analysis reveals a candidate signal of marginal 
significance with a projected orbital velocity amplitude of 
Kp = 97 km s~^, assuming a grey albedo spectrum. If gen- 
uine, this would suggest an orbital inclination close to ~ 37° , 
a planet mass Mp = 7.28 (±0.83) Mjup and a grey geometric 
albedo of p = 0.32 (±0.13), assuming Rp — 1.2 Rjup. If we 
feign complete ignorance of the value of Kp, our bootstrap 
Monte Carlo simulations give a probability ranging from 3 
to 15% that the detected feature is a consequence of spuri- 
ous noise from the analysis. When taking into account our 
prior knowledge of the system these false alarm probabilities 
drop to below 3%. 

In particular, the Class V albedo model - in which only 
the spectrum shortward of 550 nm is unaffected by Na I 
D absorption - gives a false-alarm probability of only 0.3% 
when the prior probability distribution for Kp is taken into 
account using Bayes' Theorem. However, we consider this 
is still too large an uncertainty for us to claim a bona fide 
detection. Our simulations show that a statistically unassail- 
able detection should produce a clearly visible, dark streak 
along the planet's trajectory in the trailed spectrogram, and 
no such streak is apparent even for the Class V model. 

The observations in the 2000 season were conducted 
at orbital phases optimised to produce the strongest pos- 
sible signal at a star-planet separation in velocity space 
sufficient to avoid blending problems. By adopting similar 
observing strategies on 8m-class telescopes, future reflected 
light searches of r Bootis should be able to double the effec- 
tive planetary signal contained in the WHT data described 
here, in only a small fraction of the 17 nights devoted to 
this search. Indeed, 2 optimally-phased clear nights on the 
Keckl/HIRES combination should be able to reproduce our 
results, whilst the increased efficiency of the HDS spectro- 
graph would allow Subaru to achieve very close to this depth 
of search over the same timescale. With this in mind, we 
believe that r Bootis remains a suitable target for future 
refiected light searches on 8m-class telescopes. 
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APPENDIX A: CALIBRATING THE 
MATCHED-FILTER ANALYSIS 

The purpose of incorporating a simulated planet signature 
into our analysis is two-fold. First, it allows us to ensure 
that any planetary signal, real or simulated, is maintained 
through the template subtraction and deconvolution proce- 
dures and can be recovered during the subsequent matched 
filter analysis. In doing so we can measure the degree to 
which any simulated signal is attenuated and infer that any 
real planetary signal would suffer a similar fate. Second, 
by using a suitable calibration factor, we can ensure that 
the matched filter detection for the simulated planet, ap- 
pears at the expected position in the resulting log(£o) = 
log p{Rp/a)^ vs Kp probability map, i.e. at log(eo) = —4.01 
in Fig. 4. Thus any potential detections within the real data 
(Fig. 5) would be suitably compensated for losses imposed 
by the extraction and analysis procedures. 

We model the reflected-light signal as a time sequence 
of Gaussians with appropriate velocities and relative ampli- 
tudes according to 



G{v,cl>,Kp) 



X exp 



g{(t),i) X 



Av„ 



(Al) 



where the amplitude Kp of the sinusoidal velocity variation 
is determined by the system inclination and stellar mass. 

The variable factor we use to calibrate the strength of 
the detected signal is the equivalent width (W,) of the stel- 
lar component of the composite line profile. By deconvolving 
the observed spectrum of HR 5694 with a list of the relative 
strengths of its spectral lines, we obtain a composite line pro- 
file exhibiting the broadening function that is representative 
of all the lines recorded in the spectrum. The deconvolved 
line profile is shown at Fig. Al alongside that for r Bootis 
itself. We recall that HR 5694 was chosen to best mimic the 
non rotationally broadened line profiles refiected from the 
near tidally locked planet. 

The planet signal should take the form of a faint copy 
of the stellar spectrum, as sharp as the deconvolved profile 
of HR 5694, located deep within the noise of the composite 
deconvolved residual profile of t Bootis, i.e. the deconvolved 
profile of the residual t Bootis spectrum following stellar 
subtraction. 

In effect the matched-filter analysis compares the W* 
value with the strength of the best fit Gaussian filter from 
the phased residual profiles, as at Fig. 2. Ideally, by using 
Wt = 4.45, any simulated planet signal should be recov- 
ered at the correct level within the log(eo) vs Kp probabil- 
ity map. However, analysis of the dataset following injection 
of a synthetic planet signal showed a 15% reduction in the 
equivalent width {W* = 3.85) was required to recover the 
fake grey albedo planet's signal at the correct strength. Any 
reduction in the strength of the simulated signal during the 
various processes would indeed manifest itself as a reduction 
in the log(eo) scaling factor, and would thus appear fainter 
than expected, necessitating a correction to W,. It is found 
that each of the extraction processes contributes to the sig- 
nal reduction, with the PGA fixed noise removal (Gameron 
et al. (2002) Appendix B) contributing 9% to the total loss. 



Kp km/s 




Figure Al. Deconvolved line profiles for HR 5694 and t Boo- 
tis, showing the average broadening exhibited by all lines within 
their echelle spectra. Tlie Hne profiles exhibit significantly differ- 
ent shapes due to the difference between the two stars' rotation 
rates. The normalised equivalent widths of the line profiles are 
W« = 4.45 and W« = 4.42 respectively. 



Any genuine reflected light signal should undergo a similar 
loss. In calibrating the grey albedo model (Section 8.2.1) we 
have therefore needed to correct for this signal loss during 
the extraction process. Similar, but less significant correc- 
tions (~ 10%) had to be applied to the Glass IV and Class 
V models in order to calibrate the {Rp/af values produced 
by the matched-filter analysis. 
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